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Abstract 

We address the problem of identifying the (nonstationary) quantum systems that 
admit supersymmetric dynamical invariants. In particular, we give a general expres- 
sion for the bosonic and fermionic partner Hamiltonians. Due to the supersymmetric 
nature of the dynamical invariant the solutions of the time-dependent Schrodinger 
equation for the partner Hamiltonians can be easily mapped to one another. We 
use this observation to obtain a class of exactly solvable time-dependent Schrodinger 
equations. As applications of our method, we construct classes of exactly solvable 
time-dependent generalized harmonic oscillators and spin Hamiltonians. 

Introduction 

The problem of the solution of the time-dependent Schrodinger equation, 

^j^\m) = H{tmt)). (1) 

is as old as quantum mechanics. It is well-known that this equation may be reduced to the 
time-independent Schrodinger equation, i.e., the eigenvalue equation for the Hamiltonian, 
provided that the eigenstates of the Hamiltonian are time-independent .Q The search 
for exact solutions of the eigenvalue equation for the Hamiltonian has been an ongoing 
effort for the past seven decades. A rather recent development in this direction is the 

application of the ideas of supersymmetric quantum mechanics [0. The main ingredient 
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^In this case, the adiabatic approximation is exact, 




provided by supersymmetry is that the eigenvectors of the bosonic and fermionic partner 
Hamiltonians are related by a supersymmetry transformation 0]. Therefore, one can 
construct the solutions of the eigenvalue problem for one of the partner Hamiltonians, 
if the other is exactly solvable. In general, this method cannot be used to relate the 
solutions of the time-dependent Schrodinger equation unless the partner Hamiltonians 
have time-independent eigenvectors. The aim of this article is to explore the utility of 
supersymmetry in solving time- dependent Schrodinger equation for a general class of 
time-dependent Hamiltonians. 

This problem has been considered by Bagrov and Samsonov and Cannata et al. 
for the standard Hamiltonians of the form H = p^/(2m) + V{x]t) in one dimension. 
Our method differs from those of these authors in the following way. First, we approach 
the problem from the point of view of the theory of dynamical invariants 0, 0. Dy- 
namical invariants are certain (time-dependent) operators with a complete set of eigen- 
vectors that are exact solutions of the time-dependent Schrodinger equation. We can 
easily use the ideas of supersymmetric quantum mechanics to relate the solutions of the 
time-dependent Schrodinger equation for two different Hamiltonians, if we can identify 
them with the bosonic and fermionic Hamiltonians of a (not necessarily supersymmet- 
ric) Z2-graded quantum system admitting a supersymmetric dynamical invariant. Unlike 
Refs. @] and [^, we consider general even supersymmetric dynamical invariants and use 
our recent results on the geometrically equivalent quantum systems to give a com- 
plete characterization of the time-dependent Hamiltonians that admit supersymmetric 
dynamical invariants. 

The organization of the article is as follows. In Sections 2, we present a brief review 
of the dynamical invariants and survey our recent results on identifying the Hamiltonians 
that admit a given dynamical invariant. In Section 3, we discuss the supersymmetric 
dynamical invariants. In section 4, we give a characterization of the quantum systems 
that admit a Hermitian supersymmetric dynamical invariant. In sections 5 and 6, we 
apply our general results to obtain classes of exactly solvable time-dependent generalized 
harmonic oscillators and spin systems, respectively . In section 6, we compare our method 



with that of Refs. and and present our concluding remarks. 



2 Dynamical Invariants 

By definition P, 0, a dynamical invariant is a nontrivial solution I{t) of the Liouville- 
von- Neumann equation 

p{t)=z[I{t),H{t)], (2) 

where H{t) denotes the Hamiltonian. 

Consider a Hermitian Hamiltonian H{t) admitting a Hermitian dynamical invariant 
I(t), and suppose that I(t) has a discrete spectrum^. Then, Eq. may be used to 
show that the eigenvalues of I{t) are constant and the eigenvectors |A„,a;t) yield the 
evolution operator U{t) for the Hamiltonian H{t) according to 

Uit) =J2J2<bit)\\n,a;t){\n,b;0\ . (3) 

n a=l 

Here n is a spectral label, a G {1,2, ■■ ■ ,dn} is a degeneracy label, dn is the degree of 
degeneracy of A„, the eigenvectors |A„, a; t) are assumed to form a complete orthonormal 
basis of the Hilbert space, u^hit) are the entries of the solution of the matrix Schrodinger 
equation: 

2-^«"(t) = A(tK(t), «"(0) = 1, (4) 
dt 

A(t) := S^{t)-A-{t), (5) 

and S"'{t) and A"'{t) are matrices with entries 

:= (A„, a; t\H{t)\Xn, b; t) , A^, := ^(A„, a; t| ^| A„, 6; t), (6) 

respectively, g 0]. Note that ^"(t), A"(t) are Hermitian matrices and u'^{t) is 

unitary. 

In view of Eq. , 

|^,,a;t) :=f/(t)|A„,a;0) = ^<(t)|A„,6;t) (7) 



^The generalization to a continuous spectrum is not difficult. 



are solutions of the Schrodinger equation (|T]). These solutions actually form a complete 
orthonormal set of eigenvectors of I{t). We may use this observation or alternatively 
Eq. (H) to show 

I{t) = f/(t)/(0)f/t(t) . (8) 

Now, suppose that I{t) is obtained from a parameter-dependent operator I[R] as 
I{t) = I[R{t)] where 

1. -R = {R^, R^, ■ ■ ■ , R^), R} are real parameters denoting the coordinates of points of 
a parameter manifold M; 

2. R{t) determines a smooth curve in M; 

3. /[-R] is a Hermitian operator with a discrete spectrum; 

4. (in local coordinate patches of M) the eigenvectors |A„,a; R) of /[-R], i.e., the solu- 
tions of 

/[^]|A„,a;^) = A„|A„,a;^) with a G {1, 2, ■ ■ ■ , d„} , (9) 
are smooth (single-valued) functions of R] 

5. A„ and dn are independent of R] 

6. |A„, a; R) form a complete orthonormal basis. 

In the following, we shall identify |A„, a; t) with lA^^, a; -R(t)) and express |A„, a; ^) in the 
form 

\\n,a-R) = W[R]\K,a-R{Q)), (10) 

where W[R\ is a unitary operator and = Vr[-R] defines a single- valued function of R. 
Eqs. (I) and (|10D suggest 

I[R] = W[R]IiO)W[R]^ . (11) 



For a closed curve R(t), there exists T G such that R(T) = -R(O), and the quantity 

r"(r) := Te'/d'-^-C*'^' = pe'^^" (12) 



yields the non-Abelian cyclic geometric phase |]T0|, |TT|, associated with the solution \ipn', cl', t) 



In Eq. (|T2D, T and V respectively denote the time-ordering and path-ordering operators, 
the loop integral is over the closed path R{t), and is the nondegenerate non-Abelian 



generalization of the Berry connection one- form [|1^, . The latter is defined in terms of 
its matrix elements: 

A^Ji?] ■.= t{Xn,a;R\d\Xn,b;R) , (13) 

where d = Y^- dR^ d/dR^ is the exterior derivative operator on M. If A„ is nondegenerate, 
r^{t) is just a phase factor. It coincides with the (nonadiabatic) geometric phase of 
Aharonov and Anandan [p^ . 

Next, we introduce W(t) := Vr[^(t)]. Then as discussed in Hermitian Hamiltoni- 
ans that admit the invariant 

I{t) = W{t)I{0)W{t)^ (14) 

have the form 

H(t) = W(t)Y(t)W(ty - iW(t)^W(ty , (15) 

dt 

where Y{t) is any Hermitian operator commuting with /(O). Note that according to 
Eq. (^), H{t) is related to Y{t) by a time- dependent (canonical) unitary transformation 
of the Hilbert space |]15|, |16], |], 0], namely lipit)) — » W{t)\ilj{t)) . This observation may be 



used to express the evolution operator U{t) of H{t) in the form 

Uit) = Wit)Vit) , (16) 

where V{t) := Te~'^^o y{t')dt'' jg ^i^q evolution operator for Y{t). Note that Y{t) commutes 
with /(O), therefore if /(O) has a nondegenerate spectrum, Y{t) has a constant eigenbasis. 
In this case, Y{t) with different t commute and 

r(t) =e-^/o^(0^t'. ^^7^ 



Having expressed U{t) in terms of Y{t) and W(t), we can write the solutions ([7p of the 
Schrodinger equation in the form 

\iJn.a-t) := iy(t)\/(t)|A„,a;0) = KtWW^WIAn, 0) = KlWI^n, t) , (18) 

6=1 6=1 

where V^^^{t) := (A„, a; 0|V(t)|A„, 6; 0). If Y{t) with different values of t commute, this 
equation takes the form 

|V'n,a;t) :=e-^/oS."(*')<^*'|A„,a;t), 

where := (A„, a; 0|r(t)|A„, a; 0). 

We conclude this section by emphasizing that Eqs. ([15|) and (0) are valid for any time- 
dependent unitary operator W{t) satisfying Eq. (|l^). For example, one may identify W{t) 
with the evolution operator of another Hamiltonian that admits the same invariant I{t). 
Note that in general such a choice of W{t) cannot be expressed as the image of a curve R{t) 
under a single-valued function In particular, |A„,a;t)' := iy(t)|A„; a; 0) cannot be 

written as |A„,a;_R(t)) for parameter-dependent vectors |A„,a;^) that are single-valued 
functions of R. This in turn implies that |A„,a;)f:)' cannot be used in the calculation of 
the geometric phases. 

3 Z2-Graded Systems Admitting Supersymmetric In- 
variants 



A Z2-graded quantum system |]T^ is a system whose Hilbert space is the direct sum of 



two of its nontrivial subspaces 7i±, i.e., H, = 7i+ ©7i-, and whose Hamiltonian maps H± 
to T-C±. The elements of ?i+ and ?i_ are respectively called bosonic and fermionic state 
vectors, or graded state vectors with definite grading (or chirality) and 1. Operators 
preserving the grading of the graded state vectors are called even operators. Those that 
change the grading of these state vectors are called odd operators. 

In the two-component representation of the Hilbert space, where the first component 
denotes the bosonic and the second component denotes the fermionic part of a 



state vector = + the Hamiltonian has the form 

H(.)J"*^'> " V (19) 

Here Hj^{t) : 7i+ 7i+ and H_{t) : ^ 7Y_ are Hermitian operators. They are 
respectively called the bosonic and fermionic Hamiltonians. 

Now, suppose that 7i+ = 7i_ =: 7i and consider a parameter-dependent odd operator 
Q = Q[R] and an even Hermitian operator / = I[R] that satisfy the algebra of = 1 
supersymmetric quantum mechanics [H: 

Q2 = o, [Q,/] = 0, {Q,Q^ = 2I. (20) 

In particular, suppose that in the two-component representation of the Hilbert space, 

^ 
Q = 

\d 

where d = d[R] is a linear operator. This choice of Q satisfies the superalgebra (EDf) 
provided that 

V '- 

with 



(21) 



I+:=^d^d, I.:=^ddK (22) 

As is well-known from the study of the spectral properties of supersymmetric systems, 
one can use Eqs. ([2^ ) to derive the following properties of J±. 

- /+ and /_ have nonnegative spectra with the same set of positive eigenvalues A„; 

- The degree of degeneracy d^ of A„ > as an eigenvalue of J+ is the same as its 
degree of degeneracy as an eigenvalue of J_; 

- Orthonormal eigenvectors \Xn,a,±] R) of I±[R] associated with A„ > are related 
according to 



d[R]\X^,a,+;R) = y/2Kf^vta[R]\Xn,b,-;R), (23) 

b=l 

d^[R]\K,b,-,R) = y^J2vab[R]^\Xn,a,+;R) , (24) 



a=l 



where fafe[-R] are the entries of a unitary dn x dn matrix v[E\. In particular, for a 
given orthonormal set {jA^, a, +; ^) | A„ > 0} of the eigenvectors of /+[-R], 

|A„,a,-;i?) := {2K)-^'^d[R]\K,a,+-R) 

form a complete orthonormal eigenbasis of /-[-R] for — Ker(/_[^]). Here 'Ker' 
denotes the kernel or the eigenspace with zero eigenvalue. 

Next, we introduce I{t) := I[R(t)] and I±(t) := I±[R(t)] for some curve R(t) in the 
parameter space M and demand that I{t) is a dynamical invariant for the Hamiltonian 
H{t). In view of Eqs. (D, (|19]), (P), and (|2|), is a dynamical invariant for H±{t). 

We can write /±[-R] in the form (0) by requiring d[R] to satisfy 

= W4R]d{0)W4R]\ (25) 
where := d[R{t)] and W^±[^] fulfil 

|A„,a,±;^) = iy±[^]|A„,a,±;^(0)). (26) 
In view of Eqs. (|22|) and (0), 

I±[R] = W±[R]I±{0)W±[R]K (27) 

Moreover, employing Eqs. ( p!5| ) and (|16]), we can express the Hamiltonians H±{t) and 
their evolution operators U±{t) in the form 

H±{t) = W±{t)Y±{t)W±{t)^ -iW±{t)^W±{t)\ (28) 
U±{t) = W±{t)V±{t), V±{t) ■.= Te-'^oy±(t')'it\ (29) 

where Y±{t) are Hermitian operators satisfying 

[l±(t),/±(0)]=0. (30) 

For example, we can choose Y±{t) = (1-^(0)) where Pt is a polynomial with time- 
dependent coefficients. 



In view of Eq. (ITB|) , we have the following set of orthonormal solutions of the Schrodinger 



equation for the Hamiltonians H±{t). 

dn 

|^„,a,±;t) =^K"6±W|An,&,±;t) , (31) 

b=l 

where \4"^(t) := (A„, a, ±; 0|V±(t)|A„, b, ±; 0). As we discussed above, given an eigenbasis 
|A„, a, +; t) for I+{t) we can set 

\Xn,a,—,t) = {2Xn)~^^^d{t)\Xn,a,+;t) for A„ > 0. 

This identification may be used to relate the solutions (0) according to 

dn 

a, -; t) = (2A„)-i/2 J] Kr+(^)KL W^WI^n, b, +; t) , VA„ > . (32) 

b,c=l 

For the case where Y±{t) with different t commute, Eqs. ( ^TD and (|32|) take the form 

|^n,a,±;t) = e-i'o^"±(*')'^*'|A„,a,±;t) , (33) 
|^n,a,-;t) = (2A„)-VV/o[^.\(*')-.^(0]rf*'rf(t)|^^,«,+.t), VA„ > , (34) 

respectively. Here ya±it) '■= (An, a, ±; 0|Kj-(t)|A„, a, ±; 0). 

The above construction is valid for any choice of time-dependent unitary operators 
W±{t) satisfying 

d{t) = W.{t)d{0)W+{ty . (35) 

These observations together with Eq. (H) suggests a method of generating a class of exactly 
solvable time-dependent Schrodinger equations. This is done according to the following 
prescription. 

1. Choose a Hamiltonian H^{t) whose time-dependent Schrodinger equation is exactly 
solvable, i.e., its evolution operator U+{t) is known; 



2. Choose an arbitrary constant operator do and a unitary operator W_{t) satisfying 
W^O) = 1; 



3. Set J+(0) := 4c^o/2, /-(O) := c/o4/2, and = f/+(t). Then, by construction 
I+{t) := f/_|_(t)/+(0)f/|. is a dynamical invariant for Hj^{t). It also satisfies I+{t) = 
d{t)^d{t)/2 for d{t) := W_{t)doU+{ty . Note that, in view of Eq. (^) and the 
Schrodinger equation 

tj^U+it) = H+it)U+it), (36) 
this choice of W^(t) corresponds to taking l + (t) = 0. 

4. Let Y_{t) be a Hermitian operator commuting with /-(O). Then according to 
Eq. (H), if+(t) and 

/7_(t) := W4t)Y4t)W4ty - tW4t)^ W_ity (37) 

are partner Hamiltonians, and H_{t) admits the invariant /-(t) := W^{t)I_{0)W^(ty . 

The choice W+{t) = U^{t) also implies that |A„,a, = f/_,_(t)|A„, a, +; 0) are solu- 
tions of the Schrodinger equation ([^) for H+{t). Furthermore, for all An > 0, 

\^n,a,-,t) = {2K)-'/'J2^2_it)d{t)\K,b,+-t) (38) 

b=i 

are solutions of the Schrodinger equation for the Hamiltonian H^{t). These solutions 
span H- — Ker(/_(0)) = 7Y_ — Ker((io)- Again if Y_{t) with different values of t commute, 
we have 

\^Pn,a,-;t) = (2AO-'/'e-*i'o2/a-(Orf*'rf(t)|A„,a,+;t), VA, > 0. (39) 

One can also employ an alternative construction for the Hamiltonian H_{t) in which 
one still defines I+it) according to I+(t) := U+{t)I^{0)U+{ty but uses another unitary 
operator W+(t) to express it as I+(t) = W^(t)I^{0)W-^-(ty . In this way, one may choose 
W^(t) to be the image of a curve R(t) in a parameter space M under a single- valued 
function W+ = W+[R]. This is especially convenient for addressing the geometric phase 
problem for the Hamiltonians H±{t). Following this approach, one must determine V+(t) 
according to Eq. (^), i.e., 

n(t) = W4tyH4t)W4t) - ^W^ityJ^ W^it). (40) 

One then obtains the Hamiltonian H_{t) by substituting (^) in (pSf). 



4 Partner Hamiltonians for the Unit Simple Har- 
monic Oscillator Hamiltonian 

In this section we explore the partner Hamiltonians for the Hamiltonian of the unit simple 
harmonic oscillator: 

H+ = ^{p^ + x^) =a^a + ^ = aa^ (41) 

Here p and x are respectively the momentum and position operators and a := (x + ip) / \/2. 
Let (t) be a unitary operator satisfying iy„ (0) = 1 and 

d{t) := W-{t)a^ . (42) 

Then /_|_ = Sd/2 = 0^2 is a dynamical invariant for This invariant together with 

I_(t) = ^d{t)d{ty = ^W4t)a^aW4ty (43) 

form a supersymmetric dynamical invariant. The associated 'fermionic' partner Hamilto- 
nian is given by Eq. ( P7| ) where Y_{t) is a Hermitian operator commuting with /_(0) = 
a^a/2. 

For example, let 

Y-it) = ^ (2a^a + 1) = ^ (p' + x') , (44) 

where f{t) is a real-valued function, and W^(t) = W^[9{t),ip(t)] where 

W4e,<f] := e-*^^3e-*^^2e*^^% (45) 
Ki := i(x2-p'), ir2:=-^(xp + px), i^s := ^ (x^ + p^), (46) 

6' G M, and G [0,27r). Note that Y_{t) with different values of t commute and the 
operators Ki are generators of the group SU{1,1) in its oscillator representation. The 
parameter space of the operator W- is the unit hyperboloid: 

M = {{R\ R\ R") eW" \ - {R'y - {R^f + {R'f = 1}. 



We have made the choices (|^) and (^5]) for Y_{t) and W-{t) in view of the following 
considerations. 



1. Up to a trivial addition of a multiple of identity, (^) is the most general expression 
for a second order differential operator commuting with /-(O). 

2. Every element of (the oscillator representation of the) Lie algebra of SU{1, 1) may 
be expressed as W-Y^Wl with and W- given by Eqs. (|^) and (|45|), respec- 
tively. In particular, as we show in the following, these choices lead to the most 
general expression for an invariant /_ (t) and a Hamiltonian if_ (t) belonging to (the 
oscillator representation of the) Lie algebra of 5*^7(1, 1). 

In order to compute the Hamiltonian H^{t), we substitute Eqs. (i^) and (^) in 
Eq. (^) and use the 5^(1, 1) = so(2, 1) algebra, 

[K,, K2] = -iK,, [K2, K,] = iK,, [K,, K,] = tK2, 

and the Backer- Campbell-Hausdorff formula to compute the right-hand side of the result- 
ing equation. We then find, after a rather lengthy calculation, 

3 

H.{t) = J2R\t)K,, (47) 



i=l 



where 



R\t) := smh9{t)cosip{t)[2f{t) - ip{t)]-smip{t)9{t), (48) 
R^it) := smhe{t)smip{t)[2f{t) - ip{t)] + cosip{t)9{t), (49) 
R^t) := 2cosh^(t)/(t) + [1 -cosh^(t)]^(t), (50) 

and a dot denotes a time derivative. 

As seen from Eqs. (^61) and (^Tj), H^{t) is the Hamiltonian of a time-dependent gen- 
eralized harmonic oscillator with three free functions f{t), 0{t), and ip{t). According 
to our general analysis, the corresponding Schrodinger equation is exactly solvable. The 
evolution operator is given by 

U_{t) = e-Mt)K3^-ie{t)K2^iK3Mt)-F(t)]^ 

where F{t) = f(t')dt'. Furthermore, we can use the stationary solutions of the Schrodinger 
equation for the unit simple harmonic oscillator ( ^T]) to construct solutions of the Schrodinger 



equation for H_{t). The stationary solutions for the Hamiltonian (ffTI) are 



|V^„,+;t):=e-**^-|n), (51) 

where En = n + 1/2, \n) = (?i!)^-^/^a^"|0), and |0) is the ground state vector for the unit 
simple harmonic oscillator (^if) given view of Eqs. (0), (^21), 

and (|5T|) , we have the following orthonormal solutions of the Schrodinger equation for 
H_{t). 

\^l;n,-,t) = (n + l)-i/2e-^*^"e-^^W^W_(t)at|n) = e-^^-Wg-^t^W+^Wl^^e-^^W^^ |n + 1) , 

(52) 



where Cn{t) := [t - ipit)/2]n + t/2 - 3f{t)/4:. Next, we use the identity [|TS 

e^eK,^x) = \e'/'x), 



and the expression for the propagator of the unit simple harmonic oscillator [^], namely 

f/(x,t;x',0) := {x\U{t)\x') = {27iismt)-^/^e , 

to compute the solutions ( |52| ) in the position representation. This yields 

POO 

{xl^lJn, -; t) = e-^^"W / Uix, F{t) + if{ty, x', 0)</.„+i(e^W/'a;')cia;', 



where 4>n{x) '■= {x\n) are the eigenfunctions of the unit simple harmonic oscillator Hamil- 
tonian (ETI). 



5 Partner Hamiltonians for the Dipole Interaction 
Hamiltonian of a Spinning Particle in a Constant 
Magnetic Field 

Consider the dipole interaction Hamiltonian of a spinning particle in a constant magnetic 
field: 

H = bJ3, (53) 



where b is constant (-Larmor frequency), the magnetic field is assumed to be directed 
along the z-direction, and J3 denotes the 2;-component of the angular momentum operator 
J = (Ji, J2, J3) of the particle. Let W-(t) be a unitary operator satisfying IV_(0) = 1 and 

d{t) :=iy_(t)J+, (54) 

where J± := Ji ± iJ2 = J^:- Then, in view of the identity 

[J_J+,J3] = 0, (55) 

the operator 

/+ = = J_J+/2 (56) 

is a dynamical invariant for H^. Eq. (^) follows from the su{2) = so{3) algebra, 

[Ji, J2] = iJ3, [J2,J3] = iJi, [J3, Ji] = iJ2, (57) 
satisfied by Jj, the fact that is a Casimir operator, i.e., [J^, Ji] = 0, and the relation 

J- J+ = + Jl - J3 = J2 - J3(J3 + 1) . (58) 

The invariant /+ together with 

/_(t) = \ d{t)dity = iiy_(t) J+ J_l^_(t)t (59) 

form a supersymmetric dynamical invariant. The associated 'fermionic' partner Hamilto- 
nian is given by Eq. (|37| ) where Y_{t) is a Hermitian operator commuting with /_(0) = 
J+J_/2. 

Next, we note that [J+J_, J3] = 0. This suggests that we may choose Vl(t) as a 
polynomial in J3 with time-dependent coefficients. For example, we may set 

Y4t) = f{t)Js (60) 

where f{t) is a real- valued function. With this choice of Yl, we can construct a class of 
partner Hamiltonians H^{t) for representing the dipole interaction of a spinning par- 
ticle in a time-dependent magnetic field, provided that we choose W^{t) = W-[9(t),ip(t)] 
according to 

W_[e, if] := e-^^^^g-ieJsgiv'Ja ^ ^gl^ 



where 6 G [0, vr) and ip G [0, 2ti). Note that again Y_(t) with different values oft commute, 
the parameter space of the operator is the unit sphere, and 9 and if are respectively 
the polar and azimuthal angles.^ 

The calculation of Hamiltonian H_{t) for these choices of Y_ and is similar to 
that of section 4. Substituting Eqs. (|^) and (^) in Eq. (|37[) and using the 5^(2) algebra 
([57| ) and the Backer- Campbell-Hausdorff formula, we find 

3 

H^{t) = Y,R\t)J^^ (62) 

i=l 

where 

R^{t) := sin^(t)cos^(t)[2/(t) -sin(^(t)^(t), (63) 
i?=^(t) := sin^(t)sin(^(t)[2/(t) -0(t)] +cos(^(t)^(t), (64) 
/?3(t) := 2cos^(t)/(t) + [1 -cos^(t)]^(t), (65) 

As seen from these equations, the fermionic partner Hamiltonians ( |62|) to the bosonic 
Hamiltonian (^) also belong to the Lie algebra su{2)\ they form a three-parameter fam- 
ily of dipole Hamiltonians describing spinning particles in time-dependent magnetic fields. 
The solution of the Schrodinger equation for this type of Hamiltonians has been exten- 
sively studied in the literature. A rather comprehensive list of references may be found 
in i. 

In view of Eqs. (pQ]), (|60D , and (|6TD, the evolution operator for the Hamiltonian (|62D, 
for arbitrary choices of functions /, 6', and (y?, is given by 

where F(t) := £f{t')dt'. Moreover, using the supersymmetric nature of our construction, 
we may construct a set of orthonormal solutions \iprn,~]t) of the Schrodinger equation 
for this Hamiltonian from those of the constant Hamiltonian (^). 

In order to compute these solutions, we first note that W^(t) = 1. Therefore, in view 

of Eq. (^01), l + (t) = H+(t) = bJ^. Furthermore, because J+ commutes with J3, we may 

^As described in Refs. |^], it turns out that W- as given by Eq. (|l]) fails to be single- valued at the 
south pole {6 = tt). One can alternatively change the sign of 9 on the right-hand side of ( |6l| ) , in which 
case W_ becomes single- valued for all values of and (p except for 6 — 0, i.e., the north pole. 



set 

\\m,+]t) = \j,m), (67) 

where | j, m) are the well-known orthonormal angular basis vectors satisfying 

Jz\j,m) =m\j,m), J^\j,m) = j{j + l)\j,m) , (68) 

j G {0, ±1/2, ±1, ±3/2, ■ ■ ■ } labels the total angular momentum (spin) of the particle, 
and m G {— j, — j ± 1, ■ ■ • , J — 1, j} is the magnetic quantum number. Now, in view of 
Eqs. (|56|) , (|58|) and (^), the eigenvalues of /+ are given by 



Am = {j,m\I+\j,m) = j{j + 1) -m(m± 1) . (69) 



The solutions of the Schrodinger equation for the Hamiltonian (|53| ) that are associated 
with this choice of \Xm, ±; t) are the stationary solutions 



|^m,±;t) = e-^'^+|j,m) = e-^'^J^m) . (70) 

Under the supersymmetry transformation, \ipm,+;t), with m < j, are mapped to the 
following solutions of the Schrodinger equation for the Hamiltonian (|6^). 



2 [j (j ± 1) — m(m ± 1)J 



Note that here m G {— j, — j ± 1, ■ ■ ■ , j — 2, j — 1} and we have made use Eq. ( plD and 
the relations 

= (j,"^l^+(^)|j,"^) = y-it) = (j,m|F_(t)|j,m) = /(t)m, 

J+\j,m) = \J (j - m)(j ± m ± 1) |j, m ± 1). 

Next, consider the special case of the Hamiltonians ( |62D obtained by choosing 6'=constant 
and Lf = ujt for some uj G M^, namely 

iJ_(t) = hr{t){h{t)[cos{u:t)Ji + sm{ujt)J2] + f2{t)Jz] , (72) 
r(t) := ry4/(t)2 + [a;-4/(t)]c.sin2^, 
/i(t) := sin^[2/(t)-^]/[6r(t)], ^(t) := 2 cose/(t)/[6r(t)]. 



These correspond to the dipole Hamiltonians for which the direction of the magnetic field 
precesses about the z-axis and its magnitude is an arbitrary function of time. The case of 
the constant magnitude is obtained by setting /=constant. This is the well-known case 



of a spin in a precessing magnetic field originally studied in Ref. For a more recent 
treatment see |2^, 0] . 

In Section 4, we restricted ourselves to the study of the quadratic invariants (the 
invariants that are second order differential operators). This restriction determined the 
expression for the operator Y_{t). A choice of Y_(t) which includes cubic or higher powers 
of p would lead to the fermionic partner Hamiltonians that cannot be expressed as a second 
order differential operator. An analogue of this restriction for the systems considered in 
this section is the condition that Y_{t) should belong to the Lie algebra su{2) = so{3). 
Unlike the case of the harmonic oscillators, a violation of this condition does not lead to 
any serious problem for the spin systems. For example, if we take 



Y.{t) = f{t)Js + g{t)J', 
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(73) 



but keep the same choice for W^(t), i.e., (|6ll), we are led to a class of exactly solvable 
fermionic partner Hamiltonians of the form 



H'_{t) = H_{t) + H4t) 



(74) 



where H_{t) is given by Eq. (^) and H^(t) is a general quadratic Stark Hamiltonian 



describing the quadrupole interaction of a spinning particle with the magnetic field p2 |. 
A straightforward calculation yields 

3 



git)W.it)JiW.ity = g{t)[W.it)JsW.{ty]' = g{t) 



i=l 



R\t) 

The Hamiltonian (|75|) belongs to the class of quadrupole Hamiltonians 



, (75) 



sin ^(t) cos (^(t), R^{t) := sin 9 (t) sin ^{t), R^{t) := cos 9 (t). 
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(76) 



whose algebraic and geometric structure has been studied in Refs. |2J]. In particular, 
up to trivial addition of a multiple of identity, any quadrupole Hamiltonian may be written 



in the form H{t) = J2t=o P'^^a where p" are real parameters and 

eo := J|-JV3, ei:={JiJ3 + J3Ji)/V3, 62 := (J2 + ^3^2)/V3, 
63 := {J^-Jl)/VS, 64 := (J1J2 + J2Ji)/v^. 

Furthermore, the commutators [60,6/3] =: Tq,,/? generate the group Spin{5) = Sp{2) that 
acts on the set of all quadrupole Hamiltonians, [p3| .P| 



These observations suggest that one may construct supersymmetric dynamical invari- 
ants whose bosonic and fermionic components are linear combinations of the generators 
Tq,/3, i.e., they belong to the Lie algebra of Spin{5), i.e., so(5) = sp{2). This in turn 
implies that they may be obtained from constant elements of 50(5) by 50(5) rotations. 

Next, observe that both 60 and 64 commute with J3. Therefore, in our construction of 
the partner Hamiltonians for the constant dipole Hamiltonian (|53|) , we may take Y^{t) = 
^(t)To4, where ^ is a real-valued function of time. Now, if we take 

for arbitrary functions R"^, we obtain the most general invariant I-{t) belonging to so(5). 
By construction, the corresponding fermionic partner Hamiltonians H^{t) will constitute 
a large class of time- dependent exactly solvable Hamiltonians belonging to the Lie algebra 
so(5). The explicit calculation of H_{t) requires an appropriate parameterization of the 
operator W- in terms of the coordinates of its parameter space. 



6 Discussions and Conclusion 

In this article we studied supersymmetric dynamical invariants. For a given time- dependent 
Hamiltonian H^{t), we have constructed a supersymmetric dynamical invariant I{t) and 
an associated partner Hamiltonian H-(t) such that the bosonic part of /(t) is a dynami- 
cal invariant for H^{t) and the fermionic part of I{t) is a dynamical invariant for H_{t). 

^The quadrupole Hamiltonians (^6|) also arise in the study of the adiabatic evolution of a complex 
scalar field in a Bianchi type IX background spacetime p5|] . 



We have shown how the solutions of the Schrodinger equation for H^(t) may be used to 
obtain solutions of the Schrodinger equation for H_(t). 

In order to compare our approach with those of Refs. 0, H, we note that we could 
construct an even supersymmetric invariant of the form (|2l|) by requiring the supersym- 
metric charge Q to be a dynamical invariant. It is not difficult to show that substituting 
Q in the Liouville-von-Neumann equation yields the intertwining relation 



d 



d (77) 



for the operator d. Note that this relation is only a sufficient condition for I = {Q, Q^}/2 
to be a dynamical invariant. This in turn implies that our method is more general than 
that of Refs. |[. One way to see this is to substitute Eq. (^) in Eq. ([771). Using 
Eqs. (|28|) and (^), one can then reduce Eq. (^) to 

doY+it) = Y_{t)do. 

It is not difficult to construct operators Y±{t) that commute with /±(0) but do not satisfy 
this equation. 
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